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ABSTRACT 

The solar abundance controversy inspires renewed investigations of the basic physics used to develop 
solar models. Here we examine the correction to the proton-proton reaction rate due to dynamic screening 
effects. Starting with the dynamic screening energy from the molecular-dynamics simulations of Mao 
et al., we compute a reaction-rate correction for dynamic screening. We find that, contrary to static 
screening theory, this dynamic screening does not significantly change the reaction rate from that of the 
bare Coulomb potential. 

Subject headings: equation of state - nuclear reactions, nucleosynthesis, abundances - plasmas - Sun:general 



1. INTRODUCTION 



Solar models generated w i th the lGrevesse & Noels 



d 19931) or lGrevesse & Sauvall (119981) abundances agree 
quite nicely with helioseismic inferences of the sound 
speed, the location of the base of the convection zone, 
and the helium abundance in the convection zone. 



Howev er. IXsplund et al.1 d2005l |2009t) 
d2008l 120091) . and iLudwig et alj J2009h 



Caffau et al. 



revised the 

solar abundances using three-dimensional hydrody- 
namic models of the atmosphere with improved input 
physics and non-local thermodynamic equilibrium ef- 
fects, lowering abundances by up to 1/3. Solar models 
that use the new lower abundances yield worse agree- 
ment than those that use the older abundances. This 
has led to years of heated debates over whether the 
disagreement represents a solar model problem or a 
solar abundance problem. The ongoing disagreement 
has inspired a re-examination of all aspects of solar 
models, including opacities, diffusive settling, convec- 
tive overshooting, and the possible accretion of low-Z 
material late in the Sun's evolution or mass loss early 



in the Sun's evolution (see Basu & Antia 2008; Guzik 
2008; lGuzik& Mussackll2010l for reviews of mitiga- 



tion attempts). Although these proposed adjustments 
have led to some improvement in the agreement, they 



do not satisfactorily resolve the issue. Further investi- 
gations into the basic physics of the solar interior are 
required. 

Nuclear reactions generate the energy that drives 
our Sun. Developing an accurate picture of the con- 
ditions that lead to nuclear reactions is essential in or- 
der to fully understand the inner workings of the Sun. 
With that in mind, we re-examine screening effects in 
the solar core. In this work, we focus on p-p reactions 
which produce most of the nuclear energy generated in 
the Sun. 

1.1. Nuclear reaction rates 

In this section, we derive an expression for calcu- 
lating nuclear r eaction rates following the treatment of 
Clayton! ([1968). The reaction rate per unit volume be- 
tween particles of types a and f3 is a product of the 
number densities of the particles, n a and rip, and the 
average value of the product of the relative velocity v 
times the cross section <x, 



rap 



n a np 
1 + 6 a f 



(0-v)af3 



n a na C°° 
— f- i/,(E)v(E)o-(E)dE, 

1 + O a fj Jo 



(1) 
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where if/(E) is the relative velocity distribution and 
Sap accounts for reactions of like particles. Because 
we will be dealing with the correction due to dy- 
namic screening (a ratio between the unscreened and 
screened reaction rates), we can ignore the density fac- 
tor and focus on the reaction rate per pair of particles, 



f(M,T) 



where 



3/2 



cr(E)dE, (2) 



(3) 



// is the reduced mass of the pair, and the Maxwell- 
Boltzmann distribution is used for ifr(E). The cross 
section cr(E) can be defined as a product of three sep- 
arate energy-dependent factors 



S(E) l-b 
t(E)=— exp^— 



(4) 



where b = 31.28Z (1 ,ZgA 1/2 keV 1/2 , with Z a and Z B be- 
ing the charges of the interacting ions and A is the re- 
duced atomic weight. The exponential factor in this 
expression comes from the barrier penetration proba- 
bility, the inverse energy dependence comes from the 
quantum-mechanical interaction between the two par- 
ticles, and S (E) contains the intrinsically nuclear parts 
of the probability for a nuclear reaction to occur. With 
this substitution for cr(E) , Equation|2]can be re- written 
as 



Xoo / E h \ 

S(E)exp\-—-—jdE. 



(5) 



In the non-resonant reaction case, S (E) is slowly vary- 
ing with E, so we can treat it as a constant Sq eval- 
uated at the energy where exp(-E /k B T - b/E^ 2 ) is 
maximum. Then the reaction rate per pair of particles 
(without screening) can be computed as 



A = f{n,T)S 



Jo 



exp - 



k B T a/I 



dE. 



(6) 



1.2. Electrostatic screening 



Salpeter (1954) developed a treatment to include 



the effect of static electron screening on nuclear reac- 
tion rates. Here we summarize his method which we 
will use in Section|2]as the inspiration for our calcula- 
tion of the dynamic screening correction. 



We begin by writing the total interaction energy as a 
combination of the bare Coulomb potential and a con- 
tribution from the plasma: 

ZiZ 2 e 2 



t/totalO) 



+ U(r). 



(7) 



Then consider a case in which the classical impact 
parameter r c is very small compared with the charge 
cloud radius Rq and the nuclear radius r n is much 
smaller than r c . Then the barrier penetration factor for 
r„ < r < r c depends only on the expression 



E - U(r) ■ 



Z x Z ie 2 



(8) 



For distances larger than r c , the barrier penetration fac- 
tor hardly depends on the potential. U (r) must be small 
for distances greater than Rq and approach a constant 
value Uq of the order of magnitude of Z l Z2e 2 /R D for 
small r. Then, 



Rd 



U 



«c 1, 



(9) 



where E max is the relative kinetic energy for which the 
integrand in Equation|2]reaches a sharp maximum. If 
this inequality is satisfied, U (r) can be replaced by the 
potential at the origin Uq. By examining expression|8] 
we can see that the screening potential has effectively 
increased the kinetic energy by a magnitude of Uq, so 
the cross section factors for C/ to tai = ^Coulomb + tA) are 
equivalent to the unscreened factors with energy E - 
Uq. Equation |2] can then be replaced by 

A = f(M,T)f^ Eexp(-^cr(E-U )dE. (10) 

With the change of variables E' = E — Uq and the 
approximation (E' + Uq) ~ E', the reaction rate per 
pair of particles becomes 



/DO 
E' 
U 



exp 



E' + Uq 



oiE')dE'. (11) 



Because the penetration factor for E' = -Uq is so 
small, the lower limit of the integral can be set to zero 
without significantly changing the value of the inte- 
gral. We then see that 



^screened 



exp | 



-Up 
k R T 



(12) 



illustrating that the reaction rate for the statically 
screened potential can be approximated by multi- 
plying the rate for the bare Coulomb potential by 
exp{-U /k B T). 
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Salpeter then derives Uq by solving the Poisson- 
Boltzmann equation for electrons and ions in a plasma 
under the condition of weak screening (£/totai(V) <k 
ksT). He arrives at an expression for the screen- 
ing energy that is equivalent to that of the Debye- 
Hiickel theory of dilu te solutions of electrolytes 
foebye&Huckellll923l) : 



U = - 



Z^e 2 
Rn 



(13) 



The Debye length, Rq, is the characteristic screening 
length of a plasma at temperature T with number den- 
sity n which is defined by 



R? 



e 2 (n e + riiZf) 



(14) 



For a neutral proton-electron plasma, the electron 
number density n e and ion number density n, are both 
n/2, so the Debye length is just 



Rt 



£()k B T 



(15) 



1.3. Dynamic screening 



Although Salpeter's expression accurately de- 
scribes the effect of static screening, the issue of dy- 
namic screening in the hot, dense plasma of the solar 
interior remains an open question. Dynamic screening 
occurs when the screened interaction energy of a pair 
of ions depends on the relative velocity of the pair. 
Most of the ions in the solar plasma are much slower 
than the electrons and the fastest ions. The thermal 
ions are therefore not able to rearrange themselves as 
quickly around individual faster moving ions. Since 
nuclear reactions require energies several times the av- 
erage thermal energy, the ions that are able to engage 
in nuclear reactions in the Sun are the faster mov- 
ing ions, which are not accompanied by a full static 
screening cloud. 

Salpeter's derivation uses the mean-field approach 
in which the many-body interactions are reduced to an 
average interaction that simplifies calculations. This 
technique is quite useful for calculations describing the 
average behavior of the plasma. However, dynamic ef- 
fects for the fast-moving, interacting ions in hot, dense 
plasma lead to a screened potential that deviates from 
the average value. Therefore, the mean-field approxi- 
mation is not appropriate for computing stellar nuclear 
reaction rates. Instead, we use the molecular-dynamics 



method of Shaviv & Sh aviv (1996) to model the mo- 
tion of protons and electrons in a plasma under so- 
lar conditions in order to investigate dynamic screen- 
ing in p-p reactions. The advantage of the molecular- 
dynamics method is that it does not assume a mean 
field. Nor does it assume a long-time average potential 
for the scattering of any two charges, which is neces- 
sary in the statistical way to solve Poisson's equation 
to obtain the mean potential in a plasma. 

In previous work, iMao et al. I d2009h present simu- 
lation results for the velocity-dependent screening en- 
ergy of p-p reactions in a plasma with the tempera- 
ture and density of the solar core (T = 1.6 x 10 7 K, 
p = 1.6 x 10 5 kg irT 3 ). They demonstrate that the 
static screening result does not accurately represent 
this plasma, and they compute a screening energy that 
depends on the relative kinetic energy of a pair of in- 
teracting ions (see FigureHJ. In this paper, we use their 
simulation results to compute a correction to the solar 
p-p reaction rate due to the dynamic screening they 
observe. 
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Fig. 1. — Dynamic screening energy at the turning 
point for pairs of protons with a given relative kinetic 
energy. For comparison, the static screening energy 
evaluated at the average turning point of pr oton pairs 
with e ach energy is also shown. Data from Mao et al. 
(2009). 



2. METHODS 



We begin with the calculations of lMao et al. I d2009h 
for a plasma of protons and electrons with the tem- 
perature and density of the solar core. Their Figure 
5 shows the relationship between the total interaction 
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energy at the turning point and the relative kinetic en- 
ergy of a pair of interacting protons. As in Salpeter's 
static screening derivation, we can split the total inter- 
action energy into the Coulomb and screening cloud 
contributions: 



£W0*) 



ZiZ 2 e 2 



+ U(r,v). 



(16) 



The screening energy now includes a velocity depen- 
dence, so Uo can no longer be factored out of the en- 
ergy integral as was done to obtain Equation[T2]for the 
static screening case. 

Following Salpeter's calculation for static screen- 
ing, we focus on the contribution to the interaction en- 
ergy from the screening cloud at small r. This value 
of Uo is obtained from the relationship between the 
screening energy of a pair of protons at their turning 
point and their relative kinetic energy which is shown 
in Figure Q] 

The dynamic screening energy curve is described 
by the equation 

U (E) I E \ 

-P=r =0.005-0.281 exp [-235- — , (17) 
k B T M k B T ' 



which comes from the best-fit curve for the Mao et al 
(200 9j) £gc r er,Ti(y<:, E). (Note the difference in sign from 
Mao et alJ (120091) . where the screening energy was de- 



fined as a negative contribution to the total energy, 
£to tai(r) = e 2 /r - E SCKen (r, E). In this paper, we use 
the lSalpete rl dl954l) convention, as shown in Equations 
0and[l6l) 

Now we return to Equation [TO] for the screened re- 
action rate per pair of particles. The assumptions and 
approximations used in the derivation of Equation [TO] 
for the static case will be examined and justified for 
the dynamic case in Section [3721 Replacing Uo from 
the statically screened case with Uo(E) for the dynam- 
ically screened case and using definition|4]for the cross 
section, we have 



A = f(ji, T) 
x exp 



E-Uo(E) 
E b 



S(E-U (E)) 
dE. 



k B r v(£ - mm 



(18) 



Because S (E) is a slowly varying function of energy, 
we make the approximation 



and replace this function with the constant So, as was 
done in the original reaction-rate calculation in Equa- 
tion |6] 

We can now evaluate the reaction rate per pair of 
particles using Equation|6]for the unscreened case and 



£ 



E - Uq(E) 



exp 



E 
"ksT 



b 



V(£ - Uo(E)) 



dE (20) 



for the statically and dynamically screened cases. 
Equation [T3l gives the static screening Uq and Equa- 
tion [l7]gives the dynamic screening Uo{E). Because 
all three cases contain the factor f(fi,T)So, we only 
need to compute the integrals in order to compare ra- 
tios of the two screened cases to the bare Coulomb 
potential case to obtain the correction factors for the 
p-p reaction rate. 

3. RESULTS 

Table Q] shows the results of the screening correc- 
tions for solar p-p reaction rates computed from the in- 
tegrals in Equations IB1 and l20l The statically screened 
correction shows a fairly large enhancement in the 
nuclear reaction rate. Conversely, the dynamically 
screened reaction rate is almost the same as the un- 
screened rate. 

3.1. Integrands 

How does the dynamic screening energy seen in 
Figure Q] result in a reaction rate that is so close to 
the unscreened reaction rate? To answer that question, 
we compare the components of the the integrands from 
Equations l6l and l20l Both integrands can be written in 
the general form 



F(E) = H(E)J(E), 



where 



H(E) = 



J(E) = exp 



G(E) 
E 



VG(£)' 



and 



G(E) = E-U. 



(21) 



(22) 



(23) 



(24) 



S(E-U (E))«S(E) 



(19) 



The three cases only differ in the screening energy U 
which is shown for each case in Table [TJ 
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Case 


Screening energy U 


Reaction-rate correction 


Unscreened 





1 


Statically screened 


U ° ~ Rn 


1.042 


Dynamically screened 


U Q (E) = k B T (0.005 - 0.281 exp(-2.35^)) 


0.996 



Table 1 : Screening energies and the ratio of screened to unscreened nuclear reaction rates for solar p-p reactions. 



In Figure |2(a)| we see that the dynamic G(E) ap- 
proaches the unscreened G(E) for high energies. Fig- 
ure |2(b)| shows that this leads to the dynamic H{E) 
approaching the unscreened H{E) for energies above 
~2 ksT, while the static H(E) is very different from 
both the dynamic and bare H(E). Below energies of 
~2 ksT the dynamic H(E) drops rapidly away from 
the unscreened value of 1 . However, the factor H(E) is 
multiplied by J(E) in the integrand of the reaction-rate 
equations. As seen in Figure [2(c)l J(E) is very small 
to zero below energies of ~2 kgT, damping out the re- 
gion of H(E) in which the unscreened and dynamic 
results diverge. This leads to integrands for the dy- 
namic and unscreened cases that are nearly identical, 
as seen when H(E) and J(E) are multiplied together to 
give F(E) in Figure [2(d)| The Gamow-peak-like factor 
J(E) acts as a weighting function to devalue the H(E) 
contribution of the slow pairs of ions that rarely partic- 
ipate in nuclear reactions. The faster pairs that cause 
less polarization of the surrounding plasma and there- 
fore see less screening provide the main contribution 
to the reaction-rate integral. 



Next we examine the inequality in Equation [9] 
Uo/E max <K 1. In the dynamic case, £max is the ki- 
netic energy for which the integrand of Equation l20lis 
maximum. We see from Figure [2(d)| and Uq in Figure 
[T]that the peak of the integrand F(E) occurs at energies 
for which Uo/E max <K 1. 

Finally, we address the approximation that U(r) can 
be replaced by the potential at the origin Uq- While this 
is easy to see for the static case, in the case of dynamic 
screening we do not have the screening energy in the 
form U (r, E) to examine this claim. Instead, we define 
Uo(E) to be the screening energy computed at the turn- 
ing point of the approaching protons, since this is the 
relevant U (r, E) for nuclear reaction-rate calculations. 

4. DISCUSSION OF ARGUMENTS AGAINST 
DYNAMIC SCREENING 

In light of the contentious debate over the validity 
of dynamic screening, we devote this section to a dis- 
cussion of arguments that have been made against dy- 
namic screening. 



3.2. Evaluating assumptions 

Now that we have defined the integrand F(E) for 
the dynamic case, we can return to the issue of as- 
sumptions and approximations that were justified in 
the static screening rate derivation and adopted in the 
dynamic screening rate derivation. Here we assess the 
validity of these assumptions and approximations in 
the case of dynamic screening. 

The first assumption is that r c is very small com- 
pared with the charge cloud radius Rq and that the 
nuclear radius r n is much smaller than r c , leading to 
expression [8] Although the dynamic case does not 
have a traditional static screening cloud, the inequality 
r n <k r c <k Ru is still satisfied. This can be seen in the 
definition of Ru as the distance beyond which an ap- 
preciable fraction of the nuclear charge is screened by 
the polarization charge cloud, a distance that is large 
for the dynamic case. 



4.1. Incorrect derivations 

The argument for dynamic screening has been 
damaged by several derivations of alternate screen- 
ing formulae (see, for example, Carraroet al. 1988; 



Qpher & QpherkOOOUShaviv & Shavivll996tlTsytovichl 
2000) that were subsequently s hown to be incorrect . 
Her e we discuss two examples, |C arraro e al. (1988) 
and lShaviv & Shavivld 19961) . ICarraro et alJd 1988b de 
rived a modified screening potential for fusing ions 
when the G amow velocity is greate r than the ther- 
mal velocity. Brown & Sawyerl (119971) showed that in- 
cluding processes of excitation or de-excitation of the 
plasma in an interaction with one of the fusing ions 
ex actly cancels t he dy n amic modifications prop osed 
bv lCarraro etaTl (Il988b . Ishaviv & Shavivld 1996b then 
introduced a factor of 3/2 on the screening energy. 
They arrived at this result by including the interaction 
of the the screening cloud from each fusing ion in the 
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Fig. 2. — The functions G(E), H(E), J(E), and F(E) from equations[24][22][23] andl2"Tlfor the unscreened, statically 
screened, and dynamically screened cases. 



total interactio n potential. |B mggen & Goughl i 19971) 
showed that Shav iv & Shavivl (1 1996b misinterpreted 
the thermodynamics and used an incorre ct potential in 



the Sc hrodinger equation for the system. iBahcall et al. 



J2002h summarize the problems with several differ- 
ent alternative screening formulae. However, finding 
flaws in these (and other) derivations of analytical ex- 
pressions for screening deviations does not rule out the 
effect of dynamic screening. This argument only high- 
lights the difficulty in developing a general analytical 
formalism to describe dynamic screening. 



4.2. Factorability of the distribution function is 
wrong 

For the Gibbs distribution, probabilities for mo- 
menta and coordinates are independent and cannot in- 
fluence each other. This leads to the argument that 
velocities of fusing particles cannot have an effect on 
screening because the distribution function is not fac- 
torable. However, it is clear when examining individ- 
ual ions that their relative velocity affects how close 
the ions can be to each other. This argument extends to 
ions in a plasma, where the configuration of the screen- 
ing cloud of approaching ions depends on the relative 
velocity of those ions. Over the whole system, this 
velocity-dependent effect averages out to the Gibbs 
distribution, but each screening cloud is not identical 
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to the average configuration of a screening cloud in 
that system. 

Solar nuclear reactions select a biased sample of the 
ions in the system. These nuclear reactions involve 
mainly the fastest ions, not a random sample of all ions 
in the system. Therefore the velocity distribution must 
be multiplied by the velocity-dependent screening en- 
ergy before integration instead of beginning with an 
average value of the distribution. 

4.3. Higher-order terms 

Do dynamic screening results imply that higher- 
ordered terms are required? The screening energy can 
be expressed as a power series expansion in the plasma 
coupling parameter A: 

Ell! = _^ (i _ exp (-*)) - A 2 f(x, A) , (25) 
KbL x 



where 



A 



ZiZ 2 e 2 



(26) 



R D k B T 

x = r/R D , and fix, A) is given bv lDeWittl d 1965b . The 
first term reduces to the Debye-Hiickel weak screen- 
ing result shown in Equation[T3l 

For the temperature and density of our simulations 
(T = 1.6 x 10 7 K, p = 1.6 x 10 5 kg m" 3 ), A = 0.05, 
so higher-order terms are small. Therefore, if dynamic 
screening could be described by higher-order terms in 
the expansion, the effect should be much smaller than 
the first term. However, the dynamic correction is not 
just a higher-order term in the expansion. Dynamic 
effects come from a different approach to determin- 
ing screening effects. Instead of deriving an expres- 
sion for screening based on average properties of the 
system, we examine the formation of the screening 
clouds themselves and do not average out the velocity- 
dependent nature of the clouds. 

4.4. Observational confirmation 

What observational evidence can confirm any 
screening effect in the nuclear reactions in the Sun? 
Many early discussions of dynamic screening were 
motivated by the neutrino problem which has since 
been resolved with neutrino oscillation theory. In addi- 
tion, includ ing dynamic screening co rr ections in mod- 
els wit h the Grevesse & Noels ( 1993 ) or Grevesse & Sauvall 
d 19981) abundances worsened agreement with helio- 
seismic constraints (see, for example, IWeiss et al. 



2001). However, the solar abundance problem pro- 
vides renewed motivation for exploring dynamic 
screening in solar nuclear reactions. 

Before 2005, solar models with the latest input 
physics reproduced the sound-speed profile deter- 
mined from helioseismic inversions to within 0.4% 
and also provided good agreement with the seismically 
inferred convection zone depth and convection zon e 
helium abun d ance. Then Asplund et al.l (120051, |2009), 
Caffau et alj d2008l 120091) and iLudwig et al.M 2009) 



began using three-dimensional hydrodynamic models 
of the solar atmosphere with improved input physics 
and non-local thermodynamic equilibrium effects to 
determine solar atmospheric abundances. These re- 
vised calculations lowered element abundances by up 
to 1/3. When the lower abundances are incorporated 
in solar models, the sound-speed profiles, convection 
zone depths, and convection zone helium abundances 
give worse agreement with helioseismic constraints 
than models with the old, higher abundances. Many 
attempts have been made to improve agreement by 
adjusting the physi cs or evolutionary a s sumptions in 
solar models (see iBasu & AntialEool iGuzikl 12008 : 
iGuzik & Mussackl l2010). Although these adjustments 
have shown some improvement, no model using the 
new lower abundances agrees as well with the he- 
lioseismic constraints as the models using the older 
abundances. 

In a forthcoming paper, Mussack & Guzikl ( in prep. ) 
incorporate the dynamic screening correction shown 
here for solar p-p reaction rates into solar models with 
the new lower abundances. They show that includ- 
ing this correction in solar models improves the sound 
speed discrepancy in the solar core, as shown in Figure 
|3] This improvement does not fully reconcile the new 
abundances with helioseismic constraints, but it is a 
step in the right direction. Perhaps in combination 
with other changes, dynamic screening corrections 
could contribute to a solution to the solar abundance 
problem. 

5. SUMMARY 

We have shown that dynamic screening in solar p-p 
reactions does not reproduce the enhancement of reac- 
tion rates that is predicted by Salpeter's static screen- 
in g approximation. In fact, the dynamic screening seen 
bv lMao et alj 12009) shows essentially no correction to 
the unscreened reaction rate. 

Although p-p reactions in the core are the main 
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source of nuclear energy generated in the Sun, this 
reaction-rate correction is only the beginning of un- 
derstanding how including dynamic effects will alter 
a full solar model. The reaction-rate correction must 
be generalized to treat other temperatures, densities, 
compositions, and reactions. In addition, the effect of 
dynamic screening on the equation of state must be ex- 
amined. Until we can meet both of these challenges, 
dynamic screening cannot be incorporated completely 
and consistently in solar and stellar models. 

We thank Hugh DeWitt for useful discussions and 
Dan Mao for her simulation results. This work was 
supported in part by grant AST-0708568 of the Na- 
tional Science Foundation. 
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